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Phase fluctuations of a d-wave superconducting order parameter are theoretically studied in the 
context of high-T c cuprates. We consider the t — J model describing layered compounds, where the 
Heisenberg interaction is decoupled by a d-wave order parameter in the particle-particle channel. 
Assuming first that the equilibirum state has long-range phase order, the effective action <S e // is 
derived perturbatively for small fluctuations within a path integral formalism, in the presence of 
the Coulomb and Hubbard interaction terms. In a second step, a more general derivation of <S e // 
is performed in terms of a gradient expansion which only assumes that the gradients of the order 
parameter are small whereas the value of the phase may be large. We show that in the phase-only 
approximation the resulting <S e // reduces in leading order in the field gradients to the perturbative 
one which thus allows to treat also the case without long-range phase order or vortices. Our result 
generalizes previous expressions for 5 e // to the case of interacting electrons, is explicitly gauge 
invariant, and avoids problematic singular gauge transformations. 
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I. INTRODUCTION 

The nature of the pseudogap in undcrdoped high-T c 
oxides and its relation to superconductivity remains one 
of the unsolved problems in high-T c superconductivity. 
One scenario assumes that the physics in the under- 
doped and optimally doped region is mainly determined 
by the competition of the superconducting order pa- 
rameter with a second one in the particle-hole channel. 
Possible candidates are antiferromagnetism, charge- and 
spin-density waves^*£ in particular with an internal d- 
wave symmetry*^ and stripes^ One experimental con- 
straint is that this second order parameter is strongly 
anisotropic, being large along the k x and k y axis of the 
Brillouin zone, and practically zero near the diagonal 
k x = k y . This requirement is most naturally fulfilled if 
this order parameter has e?-wave symmetry like the super- 
conducting order parameter. A second scenario assumes 
that no instabilities or strong fluctuations in the particle- 
hole channel arc relevant in the undcrdoped regime but 
that phase (and perhaps amplitude) fluctuations of the 
superconducting order parameter are important in this 
region^ 

Presently, it is unclear which of the above two sce- 
narios is more realistic in describing high-T c cuprates. 
In this paper we want to explore some aspects of the 
second approach in more detail. While it seems that 
the first scenario has become more popular than the sec- 
ond one we think that there are enough reasons to study 
phase Auctions of the superconducting order parameter. 
For instance, order parameters due to structural phase 
transitions with finite momentum transfers may easily 
be anisotropic because of strongly varying nesting prop- 
erties along the Fermi surface. However, whether one can 
achieve along this route the observed large anisotropy or 
even d-wave symmetry of the pseudogap in a generic way 
remains unclear. Experimental SIN tunneling data in the 
cuprates also show a rather symmetric one-particle den- 



sity of states, both in the superconducting and the pseu- 
dogap region, with respect to the chemical potential. 7 
Such a behavior is characteristic for superconductivity 
but not generic for densities with structural order pa- 
rameters in the particle-hole channel, especially, if they 
are associated with large momentum transfers. These 
problems vanish if one assumes that only one complex 
order parameter related to the observed d-wave super- 
conductivity is relevant. 

Fluctuations effects in superconductors arc often 
described by the time-dependent Ginzburg-Landau 
equations^ They describe a relaxational behavior of the 
order parameter back to equilibrium, i.e., their time- 
dependent part is of first order in the time derivative 
and independent of momentum. In the pure system they 
are most successfull near T c and at low frequencies com- 
pared to the temperature. In the following we want to 
consider a different regime, namely, temperatures sub- 
stantially below the mean-field transition to supercon- 
ductivity where the mean-field order parameter is large 
but varies slowly in space and time because of a small 
phase stiffness constant AiS The t-J model and its con- 
straint to have no doubly occupied sites leads necessarily 
to small values for A's at small dopings. At low temper- 
ature A is determined by the diamagnctic term which is 
given by the Fourier transform of the one-particle mo- 
mentum distribution function at non-vanishing lattice 
vectors. The constraint of having no doubly occupied 
sites implies that the maximum occupation of a momen- 
tum state is 1 + 5 instead of 2 as in the free case>i£ This 
together with the sum rule shows that the diamagnctic 
term must vanish linearly in S in the limit S — > 0. As 
a result large phase fluctuations should occur at small 
dopings. In agreement with this picture we derive in the 
present work an effective action for phase fluctuations not 
as a power but as a gradient expansion in the phase. The 
method we follow also allows to take into account inter- 
action terms between the electrons such as the Hubbard 
or the long-range Coulomb interaction. This feature is 
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important because the small phase stiffness at low dop- 
ings is caused by correlation effects. At the same time 
our treatment leads to an explicitly gauge invariant form 
for the effective action for phase fluctuations in contrast 
to many previous treatments. 



Most of the previous derivations of the effective ac- 
tion for phase fluctuations employed singular gauge 
transformation»iiii2iiii4. Such a transformation means 
in the static case for the Bogoliubov equations that 
A(r) -> A(r) • e - i<E> M, u(r) -► u(r) • e^M/ 2 , v(r) -► 



v(r) 



„i*(r)/2 



, A(r) -> A(r) - (hc/2ie)V<f>(r)J£ A,u,v 



are the pair, electron and hole wave functions, respec- 
tively, A the vector potential, and $ is an arbitrary func- 
tion. This transformation allows to remove completely 
the fluctuations in A and to describe them as fluctuations 
in A. However, A, u and v have to be unique functions 
of r, i.e., their phase can only change by multiples of 2tt 
after passing through a closed loop. This condition is in 
general violated in the presence of vortices after having 
performed such a gauge transformation. The importance 
of a possible non-uniqueness of wave functions for phys- 
ical quantities^ is presently not clear, most authors ne- 
glect this problem whereas, according to Refs. H^IT^Il^ . 
it causes very interesting effects such as non-Fermi liquid 
behavior, power laws of correlation functions, etc., in the 
normal state. In order not to violate the basic require- 
ment of uniqueness of wave functions we avoid completely 
problematic gauge transformations and derive the effec- 
tive action for phase fluctuations by means of a gradient 
expansion in the order parameter. 



The paper is organized as follows. In section II the 
Hamiltonian is specified and phase fluctuations of the 
order parameter are introduced. In section III the ef- 
fective action for phase fluctuations is derived under the 
assumption that the phases 9 deviate only little from 
a homogenous state, for instance, 6 = 0. The micro- 
scopic quantities appearing in the action arc one- and 
two-particle Green's functions associated with the Hamil- 
tonian without external potentials or phase fluctuations. 
A more general derivation for the effective action is given 
in section IV using a gradient expansion which only as- 
sumes that the gradients of the order parameter are small 
whereas 9 may be large. It is shown that the two effec- 
tive actions derived in sections III and IV are equivalent 
in the phase-only approximation for the order parame- 
ter. Section V contains a discussion of the results and 
the conclusions. In section III relations between different 
correlation functions describing density, current and pair 
fluctuations have been inferred from non-singular gauge 
transformations and used in obtaining the final form for 
the effective action. We found it useful to check these 
relations in the non-interacting case and in the long- 
wavelength, low-frequency limit directly without refer- 
ring to gauge invariance arguments. The details of this 
calculation are given in the appendix. 



II. HAMILTONIAN AND PHASE 
FLUCTUATIONS OF THE ORDER PARAMETER 

We consider a generalized t-J modelj22*2I which also 
contains repulsive Coulomb and Hubbard interactions, 
on a Bravais lattice consisting of layers of squares along 
the x and y axis. Its Hamiltonian is given by 



Hint = ^ Vij n i n 3 + fll 



(1) 

(2) 



cr is the SU(2) spin color, cj CT (ci a ) the creation (annihi- 
lation) operator of a spin a electron on the site i, J the 
Heisenberg interaction, Vij the Coulomb interaction be- 
tween the sites i and j, and U a repulsive Hubbard term. 
tij is the electronic hopping term between the sites i and 
j, {) denotes a pair index for nearest neighbor sites on 
the same layer. Si and n% are the spin and occupancy 
number operators of site i, respectively. In order to sim- 
plify the notation later we will put a = h = c = 1 , where 
a is the lattice constant of the square lattice and c the 
velocity of light. We also assume that the lattice contains 
N s sites. Let us introduce the singlet pair operators 

B h = c 1t c ]| _ 4Al ' Bi i = c J'i Ci T " c n c n » ( 3 ) 

which respectively create and annihilate a singlet on the 
bond H can then be expressed equivalently asi£ 



*j C icr C ja 



J 



(i,3) 



Hi 



(4) 



To derive the partition function, we work in a path 
integral formalism ^^24^ Heisenberg interac- 

tion term is decoupled using the Hubbard-Stratonovitch 
transformationr&Si the resulting complex fields A*, A 
correspond to the superconducting order parameter^ 
The partition function of the t — J model (0} in the imag- 
inary time formalism is given by 



Z = \ V^VVVAVA* exp( - S) 



(5) 



with the action 

r0 



S 



dr 



_ i,c 

^|i[A^(r)5*.(r)+A*.(r)^(r)] 

-S int ■ (6) 



1 

8J 



A«(r) | 2 
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Sint is the contribution to the action due to Timt , P is the 
thermal factor, (5 = /i is the chemical potential 

controlling the electron density while '5 and B* , B 
are Grassmann variables corresponding to the electronic 
operators c\ c and singlet pair operators B\ B, respec- 
tively. 

The complex bond variables A^+z and A^j+y may be 
written without loss of generality as 



Ay(r) =| A«(r) | 7 «^ M , 
-1 for j = i + x 
- 1 for j = i + y ' 



(7) 



7ij 



where a: and y are basis vectors of the direct lattice. The 
equilibrium cZ-wave order parameter value is given by 



| A y (r) |=A , & 3 -(t) = , 



(8) 



where Ao is a real number independent of i and r. In 
the present work, the fluctuations of the amplitude are 
neglected. The spatial variations of the phase <f> are sup- 
posed to be independent of the x and y directions. This 
allows us to introduce the average phase of the two bonds, 
and to put the phase variable on the lattice points, i.eH 



(9) 



where Qi is a real field. The fluctuating order parameter 
is thus assumed to be 



i[S I (r)+e, +i (r)]/2 



A m+s (t) = A • e 
A M+ . 5 (r) = -A„.c^M+^M]/2 



(10) 



It is well-known that phase fluctuations induce charge 
and current fluctuations and thus electric and magnetic 
fields^ Therefore the resulting electromagnetic field has 
to be considered in a general formalism by including 
scalar and vector potentials in the action, especially, to 
guarantee gauge invariance. In the case of lattice models, 
a minimum coupling scheme to describe the interactions 
involving the electromagnetic held is given by the Peierls 
substitution, 3 " which corresponds to 



kj ■ * 1 <t(t)*j (T (t) 



- /j, — e- Ao(n,T) 
x ^ iCT (T)* iCT (r) , (11) 

tij ■ exp ( — ie J A(r, t) ■ dl 
x # iCT (r)% CT (r) , (12) 



where Aq and A = (A x , A y ,A z ) are the scalar and vector 
potentials, respectively, and (— e) is the electron charge. 
In the case of slowly varying phase fluctuations the in- 
duced electromagnetic potentials also vary slowly in time 
and space so that wc can write^i 

J A(r,T)-dltt(r i -r j )-^{A(r j ,T)+A(r i ,T)} . 



III. EFFECTIVE ACTION FOR SMALL PHASE 
FLUCTUATIONS 

In this section we first give explicit expressions for the 
change of the action away from the d-wave saddle point 
due to slow spatial and temporal phase fluctuations and 
the corresponding induced electromagnetic field. The ef- 
fective action S e ff is defined by 



Z[A,0] 

z(°) 



expl - S e ff [A, 9} 



(13) 



using the abbreviation A = (Aq,A). Z[A,6] is the to- 
tal partition function, containing also the effects due to 
fluctuations in A and 9. 2^°' denotes the d-wave saddle 
point partition function in the presence of the Coulomb 
and Hubbard terms, and S^ ' its corresponding action 



5<°> = 



r/r 



Z<°) = J cxp( - 5 (0) ) , (14) 

^f i(7 (r){9 T -/i}*i CT (r) 

i . a 

i;{iAo7«[^(r)+fltf(r)] 

+ S int ■ (15) 



By expanding perturbatively the action 10 with re- 
spect to the phase exponential factors and electromag- 
netic potentials appearing in Eqs. (|10|) and (|ll(l . i|12[l . 
respectively, one obtains the partition function and ac- 
tion related to fluctuations in the phase and the vec- 
tor potential. We go over to imaginary frequencies and 
momentum space. The fermionic Matsubara frequencies 
arc labelled by iv m , and the bosonic ones by iuj n . It is 
convenient to define the four-dimensional wave vectors 

k = (fc Q )o<a<3 = and q = (q a )o< a <3 = (i^«,<?)- 

The Fourier transformed field variables are defined by 



* CT (kv m ,fc) = * CT (fc) 



(16) 




-if m r+ifc-T , i 



<Mt) 



It is also convenient to introduce Nambu spinor field 
operators 3 ^ 3 , by 



#(fc) = (f T (fc) *j(-fc)) , 



_ ( *t(*) 



(17) 



In this section we consider only terms up to the second 
order in the fluctuations 9. This means that we allow 
only small phase fluctuations around a homogenous state 
characterized by 9 = 0. The partition function becomes 



4 



then 



fields yields then for <S e //, denned in Eq. i|13|) . 22 



Z[A, 6] = j V^W exp( - S[A, 
S[A, 9} = 5C°) + S' A + S' e , 



(18) 
(19) 



1 

7^= E E v a {k)A a (q)$(n a , k, q) 



S'n 



^ E E rn-i,(k)A a (q)A a ,( q >) 

k.q.q 1 a.a' — l 

x $(3,fc,g + ?') , (20) 

^X>(fc,g)0(?)*(2,fc,?) 



/3N 



k,q,q' 

x l(l,fc,g+g / ) , 



with 
u (fc) = 

™aa'( fe ) 



(96 

-e , v a (k) = (~ e )-QjT for a = 1,2,3 



e 2 d 2 e fc 



2 dhj^dk^' 

^ = - E '^ iA: " 



for a, a' = 1, 2, 3 



w(fc.g) = ^ Ak + Afc +q) 



z(k,q,q') 



(A fc + 2A fc+q H 

cos(fc y )] 
$(a, ft, g) = $(fc + g)<r Q $(fc) 



Afc + q + q') 



Afc = ■^■[cos(fc : 



(21) 

(22) 

(23) 
(24) 

(25) 

(26) 

(27) 
(28) 



£fc is the hopping energy and Afc the ci-wave supercon- 
ducting order parameter. cr a denotes for a = 1,2,3 the 
Pauli matrices along the x, y, and z directions, respec- 
tively, and for a = the 2x2 identity matrix. n a is equal 
to 3 for a = and zero otherwise. The summations in 
fc-space always extend over the first Brillouin zone. 

The next step consists in performing the functional in- 
tegration over the 'J', "J fermionic fields appearing in Eq. 
(|18f) . It is achieved by noticing that Z[A, 9]/Z^ can be 
seen as a generating functional, 9, A representing exter- 
nal sources which couple to one-particle density operators 
$(a, ft, q). Abbreviating the set of variables {ai, k\, f?i}, 
{«2, &2, Q2} symbolically by 1, 2, etc., we can write 



S' A +S' e = / dlE(lMl) 



(29) 



where an explicit expression for E(l) can easily be read 
off from Eqs. I|2(J|) and l|21|) . The integration over Fermi 



S eff [A,9] = J dig c (l)E(l) 

- i / dld2 C (1; 2)E{l)E{2) + ... , (30) 

Q c are connected Green's functions. They are related to 
the usual Green's functions 0, defined by 

0(1; ...;n) = / $(l)...$(n)e- s ° /Z , (31) 



via the cumulant expansion of Q c in terms of Q obtained 
from the identitj^ 



W[A,6] =ln 



Z[A,6 



(32) 



where W[A, 9] is the generating functional for connected 
Green's functions. Explicitly, one obtains C (1) = 0(1), 
C (1;2) = 0(1; 2) - 0(1)0(2), etc. We point out that 
in the general case, these connected Green's functions 
have to be calculated in the presence of the V and U 
interaction terms. 

Expressing E in terms of A and 9 it is clear that no lin- 
ear terms in A or 9 can appear in S e ff for a non- vanishing 
momentum q. The quadratic terms in A and 9 of <S e « 
become 



Seff [A, 9) = i E A a (q)IC^(q)A a/ (-q) 
W(>m 60 {q)e{-q) 



g, a, a' 

1 
2 

J2Mq)^ e (q)8(-q) - (33) 



with 



piNs k,k' 



0N. 



x Q(n a ,k,q;n a ',k',-q) 

^™-i,(fc)(i-M(i-<W 



x 0(3, ft, 0) , 



(34) 



IC ee (q) = -^-^ W (fe,gMfc',-g) 



k,k' 



g(2,k,q;2,k',-q) 



-—Y,z(k,q,-q)g(l t k,0), (35) 



^ s fc,fc' 



x g{n a ,k,q;2,k',-q) , (36) 
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General considerations show that the Taylor expansion 
of K, ee in powers of q has no constant or linear terms, 
that of JC AS no constant terms. We therefore write 



£ 6e (q) = E fcC(«)?^ > (37) 

a, a'— 

= , (38) 



a'=0 



where the functions JC 9 ^ a , (q) and /C^, (q) approach in 
general finite values for q — > 0. Writing 



Q (g) = g a 6»(g) , 
for a = 0, 3, we obtain for S e ff 

Seff[A,6] = \ £ (^,(g)X:^(g)A,'(-«) 



(39) 



+ a (?)X:^(g)M-g) 
+ 2A a (g)/C^,(g)^(-g)) • (40) 



Eq. Q4()|l shows that iS e // actually depends not on the 
phase 9 itself but only on its gradients. The three func- 
tions IC AA , K. ee , and K, Ae are not independent from each 
other. To see this we apply a gauge transformation to 



S, 



>f.l 



, 15 



*(<?) ^<*(g) + ixa(g) , 
«(g) M<?) - 2exa(g) , 



(41) 



where both equations hold for a = 0, 1, 2, 3 and we have 
defined Xail) = QaXil)- For the following gauge trans- 
formation it is sufficient to consider both 9 and x to be 
very small. As a result we only deal with small phase 
deviations from 9 = and a possible non-uniqueness of 
wave functions cannot play any role. Invariance of S e ff 
against this transformation yields the identities 



*w(ff) = -4e a X™,( ? ) = -2iefc£,( g ) . 



(42) 



Explicit calculations for the case of non-interacting elec- 
trons are presented in the appendix and provide direct 
checks of Eq. (|42|) . Using the above relations one finally 
finds for S e ff 



q,a,a 



A a {q) + —0 a (q) 
ze 



K^(g) 



A a >(-q) + —9 a ,(-q) 
2e 



(43) 



The matrix K. AA , connects the charge and current in- 
duced by an applied external potential A. Our micro- 
scopic expression, Eq. (|34|l . just represents a generaliza- 
tion of the usual expression^ to the case with interac- 
tions between electrons. JC AA (q) is the density-density 



correlation function. For a charged system it is conve- 
nient to write it in the form, 



^ A {q) 



l + V(q)fC AA (q) 



(44) 



where JC is the irreducible part of /C, i.e., it contains 
all diagrams to K, which cannot be decomposed into two 
parts by cutting one Coulomb line, 35 - Taking into account 
the layered structure of high-T c cuprates, the Coulomb 
potential V(q) is given by^& 



V(q) 



2ire 2 d 



sinh(gii(i) 



3h(giid) — cos(q z d) 



(45) 



with gii 



1x + ^ i s the distance between layers and 



eb a background dielectric constant. IC AA approaches at 
small wave vectors the universal function 1/V(q). Sym- 
metry requires JC AA = for a! = 1,2,3. Furthermore, 
JC AA , is for a, a 1 = 1,2,3 already irreducible in the above 
sense so that no analogue to Eq. (|44|l exists in this case. 

Finally we need an equation to determine Ao. Using 
the Nambu formulation the terms in the small square 
brackett in Eq. © can be written after a Fourier trans- 
formation as 

- 7 / dT E ^ K M + K (r)] 
Jo («> 

= ^2$ k A k . k+q <S> k+q . (46) 

fc,q 

A k:k+q is, except for a factor —1/4, the Fourier trans- 
form of Ay(r)cr+ + A?.(t)<t~ with cr ± = (cri ± cr 2 )/2. 
In determining Ao we may put all phases to zero so that 
Aij(T)cr + + A*j(t)(t~ reduces to Ao7y<7i. Considering 
the deviation 8Ao from the saddle point value Ao the cor- 
responding change in the action on the right-hand side of 
Eq. H46f) yields a change linear in <SA in 5 e // after inte- 
grating out the Fermi fields as previously Requiring that 
the total linear term in <5Ao vanishes yields the desired 
equation 



Ao 
J 



■ cos 



c y )g(l,k,0). (47) 



IV. GRADIENT EXPANSION OF S, 



eff 



S e f f has been derived in section III under the assump- 
tion that the phase 9 is small. An expansion of S e ff in 
terms of gradients of the order parameter seems to be 
more satisfying because 9 is then no longer restricted to 
small values. As shown in Ref. |3f|for the case of a charge- 
density-wave state, resummations allow to transform the 
expansion of S e f t in powers of the order parameter into 
an expansion in powers of gradients of the order parame- 
ter. In the following we adapt this method to our action 
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i k+...+q ! k+... 



k+q v k+q+q 2 



k+q { +q 2 j- ^' +< ?i + <7 2 ' k+c i\ +c l 2 +q 3 



i k+...+q j ,k+...+q. +] k+...+q] 



&k+...+q. v k+...+q 




k+...+q. j, k+...+q. 



A k+...+q.,k+...+q j+l 



FIG. 1: Ring diagram of order £ in the order parameter Ak,h+ q - The solid lines denote the unperturbed Green's function Go. 



and consider first the non-interacting case Si nt = 0. The 
integration over fermions can then easily be performed 
and one obtains, dropping the constant 



A 2 



S eff {0,6] = -Trjlnt-^+A)} (48) 

1=1 

Qq is the unperturbed Green's function due to the two 
first terms in Eq. I|15|) . its expression is given in Eq. 
(|A1|) . Tr denotes the trace over k and the Nambu index. 
The £-th order term Xi in the sum over £ can be written 
in frequency-momentum space as 



9i, 



dr 



Jr-(q 1 + . ..+q t ) 



Y( qi ,q 2 ,...,q e ) , (49) 



Y(qi, <?2, qi) = ^2 Tr ' {^k,k+ qi Go(k + qi) 



x A fe+ .. 



qi-i,k+...+q t 



stands for the vector (r, r^), 
rP 



(50) 

and J dr for 

l/(/3N s ) J2i Jo dr. Tr' denotes a trace over the Nambu 
index. Xi can be visualized by a ring diagram (see Fig. 
^) where the electronic Green's function Qq (solid line) is 
scattered at the external potentials A (wavy lines). 

For the following it is convenient to write A(fc, gi) in- 
stead of Ak.k+qt ■ The first momentum k refers then to 
the relative and the second one q\ to the center-of-mass 
motion of the Cooper pair, k may assume arbitrary val- 
ues whereas gi,<72, are considered to be small. 

In the non-interacting case the calculations for 5 e // 
in section III correspond to the evaluation of ring di- 
agrams of the type shown in Fig. ^ To obtain S e jf 



in second order in 6 we had to take into account one 
and also two non-equilibrium external lines due to phase 
fluctuations yielding the second and first terms in Eq. 
(|35|) . respectively. The corresponding ring diagram of or- 
der i is obtained by evaluating the Green's functions in 
Eqs. I|34|) - (|36[) f° r non- interacting electrons and expand- 
ing them to order I — 2 or I — 1, respectively, in A(fc, 0). 
Let us first consider the first term in more detail. Its 
£-th order contribution can be illustrated by the diagram 
shown in Fig. 6A(k,qi) and 8A(k + qi,qj) are the 
two non-equilibrium lines with qi — —qj from momen- 
tum conservation. Between these lines the electrons are 



A(k,0) 



A(k,0) 



A(k,0) 




bA(k+q.,q.) 



A (*+$., 0) 



FIG. 2: Ring diagram of order I in the order parameter 
Ak,k+ q with momenta corresponding to small phase fluctu- 
ations around an equilibirum state with long-range order. 
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A^+g.,?^) 

FIG. 3: Ring diagram of order I in the order parameter Ak,k+ q 
with momenta corresponding to the leading non-local contri- 
bution in the gradient expansion. 

£ — 2 times scattered at the equilibrium order parameter 
with zero momentum. In this case the phase deviations 
from their equilibrium values must be considered as small 
so that the expansion in powers of 8 is appropriate. If 
the phase does not exhibit true long-range order but still 
varies slowly in space and time the ring diagrams of Fig. 
H have to be evaluated in a different way. In this case 
the momenta qi in the external lines have to be kept but 
one may expand the electron propagators in powers of qi. 
As shown in Refs. |3^ and |3y such an expansion gener- 
ates a gradient expansion for S e ff- In particular, it does 
not assume that 6 is small but only that the gradients of 
9 are small. Fig. [3] shows the distribution of momenta 
which yield the leading term in this expansion. Scat- 
tering at the sites i and j is accompagnied by a change 
in energy and momentum qi, whereas no change in en- 
ergy or momentum occurs at all the other sites. Keeping 
higher terms in the expansion of the electron propagator 
in terms of qi would yield contributions to S e ff which 
are at least of third order in phase gradients. 

Evaluating the diagram of Fig. [3] the sums over 
gi, ...,g,_i, g;+i, g^-i, g^+i, — qi can be immediately be 
carried out yielding products of A(fc,r), the Fourier 
transform of A(k, fci) with respect to k\. It is convenient 
to introduce a Green's function Q{k, r) by 

g-\k,r) =g^(k)-A(k,r) . (51) 

After a small rearrangement of terms one obtains for the 
diagram the expression 

X T = ^EE/^Tr'l^+^^A^g*) 

x ^-^{k + qur) 
xA(fc,-g,)}, (52) 



where Q^(k, r) is the contribution of order £ to Q. The 
sums over i, j and finally also over £ can also be performed 
yielding the following contribution to S e f / 

S eff = drTr'{g(k,r)A(k,q) 

k,q 

x £(fc + g,r)A(fc,-g)} . (53) 

Ring diagrams with only one non-equilibrium line 
<5A(fc, g) can be evaluated in a similar manner. Though 
these diagrams only involve the q = Fourier compo- 
nent of A due to energy and momentum conservation 
they yield products of phase fluctuations when passing 
from order parameter to phase fluctuations. Evaluating 
these ring diagrams in form of a gradient expansion yields 
the following contribution to S e f / 

S eff = E f dr Tr'{e(fc,r)A(fc,0)} • (54) 
k 

Eqs. (|5^|l and (|5*lj) represent time and space averages of 
the corresponding homogenous action where the Green's 
functions contain the local instead of the global gap. This 
gives a simple recipe to generalize an expression for S e f / 
derived under the assumption of small phase fluctuations 
to one which is valid also for large phase fluctuations in 
the leading order of a gradient expansion: one writes 
S e ff as a density in space and time and then uses at a 
given point r in space and time the homogenous expres- 
sion for S e ff with the constant gap value A(k,r) in the 
Green's functions. In the phase-only approximation 
order parameter fluctuations are solely due to fluctua- 
tions in the phase. Because the Green's functions in Eq. 
(I53|) refer for a given r to a constant phase this phase can 
be gauged away by a global gauge transformation with- 
out any additional contribution to the vector potential 
A. The right-hand side of Eq. I|53|> is clearly invariant 
against a global gauge transformation. This is true even 
separatly for the product of the order parameters and 
the susceptibility because the latter involves the same 
number of creation and annihilation operators so that 
global phases cancel. In the phase-only approximation 
we thus may put the phase in Eq. I|53|) to zero. This 
means that the ring diagram of Fig. 21 reduces to a dia- 
gram of second-order in A where the Green's functions 
(solid lines) are to be calculated with the equilibrium or- 
der parameter Eqs. (J7J) and (|SJ). This result proves that 
within the phase-only approximation the effective action 
derived for small fluctuations in section III is in leading 
order in field gradients identical with that of the gradient 
expansion. In the general case where also the amplitude 
Ao varies in time and space Eq. l|53|l clearly differs from 
the corresponding lowest-order expression in section III. 

It is evident that the above results also hold for inter- 
acting electrons. For a given skeleton diagram in inter- 
action and electron lines there is again a one-to-one cor- 
respondence between diagrams in ^-th order in A for or- 
dered and disordered ground states. For ordered ground 
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FIG. 4: (Color online) Small and large fluctuations in the 8-t 
plane. 

states this diagram is obtained by expanding electron 
lines in A(fc, 0) so that the total order is £ — 2. In the dis- 
ordered case the same diagram is obtained by expanding 
S e f f in Eq. (|3U[) up to the order £, and associating in all 
possible ways two vertices with and the remaining £ — 2 
without energy and momentum changes. Writing the en- 
ergy and momentum conservation as a Fourier integral 
over r, the momentum integration at the I — 2 vertices 
without energy and momentum changes can be carried 
out yielding products of A(fc, r). As a result a strict one- 
to-one correspondence between diagrams of the ordered 
and disordered cases is established from which the above 
recipe follows. 



ferential equation which may have besides of trivial con- 
stant also vertex solutions illustrated by the solid line in 
Fig.H 

S e ff reduces in the static limit to the Ginzburg- 
Landau form for the phase-dependent part of the free 
energy of a superconductor. For non-interacting elec- 
trons and A = it also coincides with the expression 
given in Ref. ^3- In contrast to the latter reference the 
expression of S e tf we have obtained, displayed in Eq. 
(|4*51) . is explicitly gauge- invariant due to the inclusion of 
A. Another reason why A should be kept in S e f / is that 
current fluctuations produce a magnetic field which can- 
not be gauged away. Heuristically, the above expression 
for S e ff can be obtained in the following simple way: 
Discard order parameter fluctuations in the original mi- 
croscopic Lagrangian but keep the potential A. Inte- 
grate then over the Fermi fields and make the resulting 
S e ff gauge-invariant by applying the gauge transforma- 
tion, Eq. (|41|l . and interprete the field — 2e\ as the phase 
field of the order parameter. Finally, we stress that our 
derivation of S e f / did not make use of any singular gauge 
transformation as in Refs.lT^ andll^. According to Refs. 
mum the use of such a transformation is problematic 
because it produces in general additional cut contribu- 
tions. The fact that our S e ff agrees for A = and non- 
interacting electrons with that of Ref. 0, however, sug- 
gests that such cut contributions must cancel out in S e f / 
and that a simple use of singular gauge transformations 
leads to the correct result for S e f / • 
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V. DISCUSSION AND CONCLUSIONS 

Eq. g3 together with Eqs. |33I)-(EHI) represent a 
microscopic expression for the effective action of phase 
fluctuations in a d-wave superconductor with interacting 
electrons. In section III this expression was derived under 
the assumption that 9 is small. The dotted line in Fig. 
0] illustrates this case for the T-dependence of 9. Within 
the interval [0, (3] 9 performs small oscillations around the 
equilibrium value 9 = 0. Similar pictures can be drawn 
for paths in r = (r, r) space. The gradient expansion 
in section IV allows to consider more general paths of 
the form 9(r) = 9o(r) + 89(r). 9o(r) is assumed to be 
slowly varying with r whereas 59(r) must be small, as 
illustrated by the solid and dashed lines, respectively, in 
Fig. El Since in the phase-only approximation the gradi- 
ent expansion leads to the same expression for S e f / as the 
perturbation theory Eq. l|4"5|) is actually valid for all paths 
shown in Fig. To determine S e ff for the solid curve 
it is sufficient to take in Eq. Ij43[) the long- wavelength, 
low- frequency limit in K AA ,(q). The important slowly 
varying extremal paths satisfy then a second-order dif- 



APPENDIX A: INTERRELATIONS BETWEEN 
DIFFERENT KERNELS K. IN S eff FOR 
NON-INTERACTING ELECTRONS 



In this appendix we check the relation between the ker- 
nels K AA and K given by Eq. (|42|l for an isolated layer 
and non-interacting electrons. The free Green's function 
matrix is then given by 



Go(k) 



1 



K) 2 + (&) 2 + (Afc) 2 
—Afc w m — £fc 



with the following energies 
Afc 



^-[cos(k x ) - cos(fc. y )] , 



(Al) 



(A2) 



£fc = — 2t[cos(k x ) + cos(fc y )] — At [cos(k x ). cos(fc y )] , 

where e/- includes nearest and next-nearest neighbor hop- 
ping contributions. Eq. I|47|) yields the BCS gap equa- 
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y o:Q 34 i 37 w here E k is the quasiparticle energy, E k = 

2 y/ {£,k) 2 + (A&) 2 , and the summation in fc-space is 

1_ 1 y [cos(fc a ) - cos(fc y )] tanh ( P E k \ (A3 , extended over the first (square) Brillouin zone. Per- 

J N s ^— J 2E k \ 2 / ' forming the summation over the fermionic Matsubara 



frequencies we find for the kernel K of Eq. (|35 



k \ 



{l - n F r>(E k - q ) - n FD (E k )} £ k £k-q + A fc A fc _ q 

1 



+ 



lU! n — E k — E k - q \LO n + E k + Ek-q 

{n FD (E k _ q ) - n FU (E k )} f _ £, k £k-g + A fc A fc _q ) 

2 [ EkEk-q J 



+ 



1 



1 



lUJ n — Ek + Ek-q ILOn + -Efc _ Ek-q 



(A4) 



with npD the Fermi-Dirac distribution function. 

To perform a first check of Eq. H42f> . we consider the 
frequency dependence of the kernels K, 98 and JC AA . With 
Eqs. and (|X4|) we have 



/C^(ko„,0) 



i-£{l_2n FD (i? fc )} 
(Afe) 2 



(A5) 



£ fc {(i^„) 2 -4(£ fc ) 2 } 

Using Eq. (|34f> we get for the electromagnetic field kernel 

:,2 



/C A 4 (i^,0) = ^-^Tr'[^o(fc)<T3^(fc)] , (A6) 



which gives after summation over the fermionic frequen- 
cies 



K&Vn.O) = ^]T{l-2n FD (£ fe )} 

(Afc) 2 



(A7) 



Sfe{K l ) 2 -4(S fc )2} 

By comparing Eqs (|A5() and (|A7(I one can see immcdiatly 

/C AA (iu;„, 0) = -4e 2 /C^(io; n , 0) , (A8) 

in agreement with Eq. (|42|1 . for any frequency iw n . 

The low-frequency, long-wavelength limit of the effec- 
tive action 5 e // [0, 9] can be investigated by expanding 
quadratically its kernel. We start by considering the zero- 
temperature case. The low-frequency expansion up to the 



second-order in iw„ of K. , Eq. (|A4|) . is straightforward, 
we have 



^Vn-O ig = 0) = -^£^(i Wn ) 2 . (A9) 



However, the low-momentum expansion up to the second- 
order in the q components requires lengthy algebra. 
From Eq. (|A4|) we get 



)C ee (iu n = 0,q^O) 



(A10) 



k a=x,y 



a(A fc ) 2 


^ 2 a\ 


fe) 2 A fc 


^ 2 A fc \ 


(Ek) 3 






V ^fc 2 ) 



2(A fc ) 2 - 



+ (a) 

— 2£ k A k 



2/ ^Afc 



dk a 



The previous expression can be simplified by making use 
of the BCS equation, Eq. (|A3|) . and its second derivatives 
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with respect to k a . This yields 



followed^ 



IC se (iu n = 0,q^ 0) 



(All) 



1 



4N S ^ ^ 

k a=x,y 



(la) 





'd 2 A fc \ 


, (A fc ) 2 


%\ 


E k 1 


v dkl ) 


(^) 3 


\dk a J 



(q) 2 /gAfc 

(£ fc ) 3 ^ afc« 



- 2 



(Ek) 3 \dk a J\dk a 



We know that in the case of an infinite two-dimensional 
square lattice, the summation over k in the first Brillouin 
zone can be computed by performing the replacement 



i 



N s ^ (2tt)2 

therefore we have the identity 
J_ fd 2 E k 



dkl 



= 



(A12) 



By expressing the second order derivative of E k in terms 
of and A*. , one can use the previous equation to sim- 
plify once more K." , Eq. I|A11|) . We obtain 



(A13) 



4N S ^ ^ E k \ dk 2 a 

h (X—x,y x Lt 



By combining Eqs. (|A9|) and (|A13|I . the low- frequency, 
long-wavelength limit of JC ee can finally be expressed as 



£ flfl (g->0) (A14) 
i|xo(iw„ -» 0,q = 0) • (iw n ) 2 

+ ^ A^i^O^O)-^) 2 ! • 



A similar procedure can be repeated in the finite tem- 
perature case. The low- frequency, long-wavelength limit 
of JC ee is given again by Eq. I|A14|) . with the mean field 
density-density correlation function xo and mean field 
phase stiffness Ag, 



Xo(iw n -> 0,q = 0) 



(A15) 



N 



A^iw,, = 0,q 
1 

~ N 



(A16) 



fe 



|l-2n FD (£; fc )}- 



£ fc V 9fc 2 



+ 2(3 ■ n FD (E k ){l- n Fr> {E k )} 

Our expressions 1A15II - 1 A 1611 h ave been compared 
with Eqs. (16) - (19) of Ref. £1 taken in the low- 
frequency, long-wavelength limit we consider. We obtain 
the same results. Therefore we can conclude that K, AA 
and K, are equal, up to a constant factor proportional 
to the electronic charge. It provides an additional non- 
trivial check of Eq. 
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